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ABSTRACT 


This short technical paper advocates a bootstrapping algorithm from 
which we can form a Statistically reliable opinion based on limited clinically 
observed data, regarding whether an osteo-hyperplasia could actually be 
a case of Ewing’s osteosarcoma. The basic premise underlying our 
methodology is that a primary bone tumour, if it is indeed Ewing’s 
osteosarcoma, cannot increase in volume beyond some critical limit 
without showing metastasis. We propose a Statistical method to 
extrapolate such critical limit to primary tumour volume. Our model does 
not involve any physiological variables but rather is entirely based on time 
series observations of increase in primary tumour volume from the point 
of initial detection to the actual detection of metastases. 
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INTRODUCTION 


To date, oncogenetic studies of EWS/FLI-11 induced malignant transformation 
have largely relied upon experimental manipulation of Ewing’s bone tumour cell 
lines and fibroblasts that have been induced to express the oncogene. It has been 
shown that the biology of Ewing’s tumour cells in vitro is dramatically different 
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between cells grown as mono-layers and cells grown as anchorage-independent, 
multi-cellular spheroids (MCS). The latter is more representative of primary 
Ewing's tumour in vivo (Lawlor et. al, 2002). MCS are clusters of cancer cells, used 
in the laboratory to study the early stages of avascular tumour growth. Mature 
MCS possess a well-defined structure, comprising a central core of necrotic i.e. 
dead cells, surrounded by a layer of non-proliferating, quiescent cells, with 
proliferating cells restricted to the outer, nutrient-rich layer of the tumour. As such, 
they are often used to assess the efficacy of new anti-cancer drugs and treatment 
therapies. The majority of mathematical models focus on the growth of MCS or 
avascular tumour growth. Most recent works have focused on the evolution of 
MCS growing in response to a single, externally-supplied nutrient, such as oxygen 
or glucose, and usually two growth inhibitors (Marusic et. al., 1994). 


Mathematical models of MCS growth typically consist of an ordinary differential 
equation (ODE) coupled to one or more reaction-diffusion equations (RDEs). The 
ODE is derived from mass conservation and describes the evolution of the outer 
tumour boundary, whereas the RDEs describe the distribution within the tumour of 
vital nutrients such as oxygen and glucose and growth inhibitors (Dorman and 
Deutsch, 2002). However studies of this type, no matter how mathematically 
refined, often fall short of direct clinical applicability because of rather rigorous 
restrictions imposed on the boundary conditions. Moreover, these models focus 
more on the structural evolution of a tumour that is already positively classified as 
cancerous rather than on the clinically pertinent question of whether an initially 
benign growth can at a subsequent stage become invasive and show metastases 
(De Vita et. al., 2001). 


What we therefore aim to devise in our present paper is a bootstrapping algorithm 
from which we can form an educated opinion based on clinically observed data, 
regarding whether a bone growth initially diagnosed as benign can subsequently 
prove to be malignant (i.e. specifically, a case of Ewing’s osteosarcoma) . The 
strength of our proposed algorithm lies mainly in its computational simplicity — our 
model does not involve any physiological variables but is entirely based on time 
series observations of progression in tumour volume from the first observation 
point till detection of metastases. 


LITERATURE SUPPORT 


In a clinical study conducted by Hense et. al. (1999), restricted to patients with 
suspected Ewing’s sarcoma, tumour volumes of more than 100 ml and the 
presence of primary metastases were identified as determinants of poor prognosis 
in patients with such tumours. Diagnoses of primary tumours were ascertained 
exclusively by biopsies. The diagnosis of primary metastases was based on 
thoracic computed tomography or on whole body bone scans. It was observed that 
of 559 of the patients (approx. 68% in a total sample size of 821) had a volume 
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above 100 ml with smaller tumours being more common in childhood than in late 
adolescence and early adulthood. Extensive volumes were observed in almost 
90% of the tumours located in femur and pelvis while they were less common in 
other sites (p < 0.001). On average, 26% of all patients were detected with 
clinically apparent primary metastases. 


The detection rate of metastases was markedly higher in patients diagnosed after 
1991 (p < 0.001). Primary metastases were also significantly more common for 
tumours originating in the pelvis and for other tumours in the Ewing's family of 
tumours (EFT); mainly the peripheral neuro-ectodermal tumours (PNET); (p < 
0.01). Tumours greater than 100 ml were positively associated with metastatic 
disease (p < 0.001). Multivariate analyses, which included simultaneously all 
univariate predictors in a logistic regression model, indicated the observed 
associations were mostly unconfounded. Further it has been found that the 
metastatic potential of human tumours is encoded in the bulk of a primary tumour, 
thus challenging the notion that metastases arise from sparse cells within a 
primary tumour that have the ability to metastasize (Sridhar Ramaswamy et. al., 
2003). These studies lend credence to our fundamental premise about a critical 
primary tumour volume being used as a classification factor to distinguish between 
benign and potentially malignant bone growth. 


STATISTICAL MODELLING METHODOLOGY 


Assuming that the temporal drift process governing the progression in size of a 
primary Ewing tumour of the bone to be linear (the computationally simplest 
process), we suggest a straightforward computational technique to generate a 
large family of possible tumour propagation paths based on clinically observed 
growth patterns under laboratory conditions. In case the governing process is 
decidedly non-linear, then our proposed scheme would not be applicable and in 
such a case one will have to rely on a completely non-parametric classification 
technique like e.g. an Artificial Neural Network (ANN). 


Our proposed approach is a bootstrapping one, whereby a linear autoregression 
model is fitted through the origin to the observation data in the first stage. If one or 
more beta coefficients are found to be significant at least at a 95% level for the 
fitted model then, in the second stage, the autoregression equation is formulated 
and solved as a linear difference equation to extract the governing equation. 


In the final stage, the governing equation obtained as above is plotted, for different 
values of the constant coefficients, as a family of possible temporal progression 
curves generated to explain the propagation property of that particular strain of 
tumour. The critical volume of the primary growth can thereafter be visually 
extrapolated from the observed cluster of points where the generated family of 
primary tumour progression curves shows a definite uptrend vis-a-vis the actual 
progression curve. 
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If no beta coefficient is found to be significant in the first stage, a non-linear 
temporal progression process is strongly suspected and the algorithm terminates 
without proceeding onto the subsequent stages, thereby implicitly recommending 
the problem to a non-parametric classification model. 


The mathematical structure of our proposed model may be given as follows: 


Progression in primary Ewing tumour size over time expressed as an n-step 
general autoregressive process through the origin: 


Se =o ty eS ee ee | 
. E i (1) 


Formulated as a linear, difference equation we can write: 

Det Pat ps tee Oe, Al) 

Taking S; common and applying the negative shift operator throughout, we get: 
[-1+ BE'+B,E°+..+B,E "|S, =-€ ain 

Now applying the positive shift operator throughout we get: 


|-E" + p E” 7 p 2 Toot p,\S, = (IV) 


The characteristic equation of the above form is then obtained as follows: 


-r+ pr + pa"? tat, +E)=0 w 


Here r is the root of the characteristic equation. After solving for r, the governing 
equation can be derived in accordance with the well-known analytical solution 
techniques for ordinary linear difference equations (Kelly and Peterson, 2000). 


SIMULATED CLINICAL STUDY 


We set up a simulated clinical study applying our modeling methodology with the 
following hypothetical primary Ewing tumour progression data adapted from the 
clinical study of Hense et. al. (1999) as given in Table | below: 
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Table I 
Observation (t) Primary Ewing tumour volume (in ml.) 


(At point of first detection) 


(At the point of detection of metastasis) 





Figure | 


Primary Tumor Growth Curve 
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The temporal progression path of the primary growth from the point of first 
detection to the onset of metastasis is plotted above in Figure I. 


We have fitted an AR (2) model to the primary tumour growth data as follows: 
E (S,) = -1.01081081S,., + 5.32365561S,., (VI) 


The R? of the fitted model is approximately 0.8311 and the F-statistic is 9.83832 
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with an associated p-value of approximately 0.04812. Therefore the fitted model 
definitely has an overall predictive utility at the 5% level of significance. The 
residuals of the above AR (2) fitted model are given in Table II as follows: 


Table Il 


Observation Predicted S Residuals 


-5,05405405 12.05405405 


19.5426024 -10.5426024 


28.168292 -9.168292003 


28.7074951 10.29250488 


61.7278351 29.27216495 


| 15.638785 -13.63878518 





The average of the residuals comes to 3.044841. Therefore the linear difference 
equation to be solved in this case is as follows: 


X = -1.01081081X,., + 5.32365561X;,.. + 3.044841 (VII) 


Applying usual solution techniques, the general solution to equation (VII) is 
obtained as follows: 


X: = c4 (2.43124756) * + cp (-3.44205837) | (VIII) 


Here c, and cz are the constant coefficients which may now be suitably varied to 
generate family of possible primary tumour progression curves as in Figure Il 
below: 


Figure I] 
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In the above plot, we have varied cz in the range 0.01 to 0.10 and imposed the 
condition c4 = 1 — c2. The other obvious condition is that choice of c4 and c> would 
be such as to rule out any absurd case of negative volume. Of course the choice of 
the governing equation parameters would also depend on specific clinical 
considerations (King, 2000). 


CONCLUSION 


From Figure Il, it becomes visually apparent that continuing increase in the 
observed size of the primary growth beyond approximately 52 ml. in volume would 
be potentially malignant as this would imply that the tumour would possibly keep 
exhibiting uncontrolled progression till it shows metastasis. This could also be 
obtained arithmetically as the average volume for t = 5. Therefore the critical 
volume could be fixed around 2 ml. as per the computational results obtained in 
our illustrative example. 


Though our computational study is intended to be purely illustrative as we have 
worked with hypothetical figures and hence cannot yield any clinical conclusion, 
we believe we have hereby aptly demonstrated the essential algorithm of our 
statistical approach and justified its practical usability under laboratory settings. 
This is not intended to be a clinical paper and in defense, it may be stated that 
given the immense potential of computer modeling and simulation in present times, 
a number of mathematical tumor growth studies are being conducted now that are 
essentially theoretical rather than clinical in nature (Swanson et. al. 2003; 
Tabatabai, Williams and Bursac, 2005). 


We have used a difference equation model rather than a differential equation one 
because under practical laboratory settings, observations cannot be made 
continuously but only at discrete time intervals. However, as already stated 
previously, the numerical example we have provided is purely demonstrative - its 
purpose is illustrating the computational modeling technique rather than drawing 
clinically viable conclusions based on the model coefficients. Having said that, we 
definitely feel that there is an immediate scope of taking our line of research further 
forward by actually implementing an autoregressive process to model in vitro 
growth of MCS with clinically observed data. 
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